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Shattering of dust grains in the interstellar medium is a viable mechanism of small grain production in galaxies. 
We examine the robustness or uncertainty in the theoretical predictions of shattering. We identify Pi (the critical 
pressure above which the deformation destroys the original lattice structures) as the most important quantity in 
determining the timescale of small grain production, and confirm that the same Pi jt (t is the duration of shattering) 
gives the same grain size distribution [n(a), where a is the grain radius] after shattering within a factor of 3. 
The uncertainty in the fraction of shocked material that is eventually ejected as fragments causes uncertainties in 
n(a) by a factor of 1.3 and 1.6 for silicate and carbonaceous dust, respectively. The size distribution of shattered 
fragments have minor effects as long as ctf < 3.5 (the size distribution of shattered fragments oc a _Qf ), since the 
slope of grain size distribution n(a) continuously change by shattering and becomes consistent with n(a) oc a~ 3 5 . 
The grain velocities as a function of grain radius can have an imprint in the grain size distribution especially for 
carbonaceous dust. We also show that the formulation of shattering can be simplified without losing sufficient 
precision. 

Key words: cosmic dust, galaxy evolution, grain size distribution, interstellar medium. 



1. Introduction 

The evolution of dust in galaxies is important to the under- 
standing of galaxy evolution, since dust grains govern some 
fundamental physical processes in the interstellar medium 
(ISM). First, they dominate the absorption and scattering of 
the stellar light, affecting the radiative transfer in the ISM. 
Second, the dust surface is the main site for the formation 
of molecular hydrogen. The former process is governed by 
the extinction curve (absorption and scattering coefficient as 
a function of wavelength; Hoyle and Wickramasinghe, 1969; 
Draine, 2003) and the latter by the total surface area of dust 
grains (Yamasawa et al, 2011). Since the extinction curve 
and the total grain surface area both depend strongly on the 
grain size distribution, clarifying the regulating mechanism 
of grain size distribution is of particular importance in un- 
derstanding those important roles of dust. 

Mathis et al. (1977, hereafter MRN) show that a mixture 
of silicate and graphite with a grain size distribution [num- 
ber density of grains per grain radius, denoted as n(a) in this 
paper] proportional to a~ 35 (a grain size distribution with a 
power index of —3.5 is called MRN grain size distribution), 
where a is the grain radius (a ~ 0.001-0.25 /im), repro- 
duces the Milky Way extinction curve. Pei (1992) shows 
that the extinction curves in the Magellanic Clouds are also 
explained by the MRN grain size distribution with different 
abundance ratios between silicate and graphite. Kim et al. 
(1994) and Weingartner and Draine (2001) have applied a 
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more detailed fit to the Milky Way extinction curve in or- 
der to obtain the grain size distribution. Although their grain 
size distributions deviate from the MRN size distribution, the 
overall trend from small to large grain sizes roughly follows a 
power law with an index near to —3.5. In any models that fit 
the Milky Way extinction curve, the existence of a large num- 
ber of small (a < 0.01 ^m) grains is required. The existence 
of such small grains is further supported by the mid-infrared 
excess of the spectral energy distribution in the Milky Way 
(e.g., Desert et al, 1990; Draine and Li, 2001). 

Dust grains formed by stellar sources, mainly supernovae 
(SNe) and asymptotic giant branch (AGB) stars, are not suf- 
ficient to explain the total dust mass in the Milky Way ISM, 
because the timescale of dust enrichment by these sources is 
much longer than that of dust destruction by SN shocks (e.g., 
McKee, 1989; Draine, 1995). Moreover, both SNe and AGB 
stars are suggested to supply large (a > 0.1 /im) grains into 
the ISM and cannot be the dominant source of small grains 
(Nozawa, et al, 2007; Hofner, 2008; Mattsson and Hofner, 
2011; Norris et al, 2012). Thus, some interstellar processes 
(or more precisely, non-stellar processes) are necessary to 
explain the large abundance of small grains. 

There are some possible processes that efficiently mod- 
ify the grain size distribution. Hellyer (1970) shows that 
the collisional fragmentation of dust grains finally leads to 
a power-law grain size distribution similar to the MRN size 
distribution (see also Bishop and Searle, 1983; Tanaka et al, 
1996). Hirashita and Yan (2009) show that such a fragmen- 
tation and disruption process (or shattering) can be driven 
efficiently by turbulence in the diffuse ISM. Dust grains are 
also processed by other mechanisms. Various authors show 
that the increase of the total dust mass in the Milky Way ISM 
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is mainly governed by grain growth through the accretion of 
gas phase metals onto the grains (we call elements compos- 
ing dust grains "metals") (e.g., Dwek, 1998; Zhukovska et 
al., 2008; Inoue, 2011; Asano et al, 2012). Grain growth 
mainly works in the dense ISM, such as molecular clouds. 
In the dense ISM, coagulation also occurs, making the grain 
sizes larger (e.g., Hirashita and Yan, 2009; Ormel et al., 
2009). In the diffuse ISM phase, interstellar shocks associ- 
ated with supernova (SN) remnants destroy dust grains, espe- 
cially small ones, by sputtering (e.g., McKee, 1989; Nozawa 
et al., 2006). Shattering also occurs in SN shocks (Jones 
et al., 1996). All these processes above other than shat- 
tering do not contribute efficiently to the increase of small 
(a < 0.01 /im) grain abundance. Indeed, Asano et al. (2013) 
show that, unless we consider shattering, the grain size dis- 
tribution is biased to large (a ~ 0.1 /im) sizes. Therefore, 
shattering is important in reprocessing large grains into small 
grains. 

Based on the formulation developed by Hirashita and Yan 
(2009), Hirashita et al. (2010, hereafter H10) consider small 
grain production by shattering in the context of galaxy evo- 
lution. H10 assume that Type II SNe (SNe II) are the source 
of the first dust grains, which are biased to large sizes (> 
0.1 /im) because small grains are destroyed in the shocked 
region in the SNe II (Nozawa et al., 2007, hereafter N07). 
Starting with the size distribution of SNe II dust grains in 
N07, H10 solve the shattering equation to calculate the evo- 
lution of grain size distribution. They show that the velocity 
dispersions acquired by the dust grains in a warm ionized 
medium (WIM) is large enough for shattering to produce 
a large abundance of small grains on a short (< 10 Myr) 
timescale. As mentioned above, even if we consider other 
sources of dust such as the dust formation in the wind of 
AGB stars and the accretion of metals onto grains in the in- 
terstellar clouds, shattering is generally necessary to produce 
small grains (Asano et al. 2013). 

Considering that shattering is a unique mechanism that 
produces small grains efficiently in the ISM, it is important 
to clarify how robust or uncertain the theoretical calculations 
of shattering are. Since there are some basic physical param- 
eters regulating shattering, it is crucial to clarify how they 
affect the grain size distribution or to which parameter shat- 
tering is sensitive. Although any models should commonly 
contain those basic parameters, there could be some uncer- 
tainties in the formulation itself, or there could be unneces- 
sary complexity which cannot be constrained from observa- 
tions anyway. Thus, in this paper, we examine the robustness 
of shattering calculations for the interstellar dust by chang- 
ing some major parameters in a model and comparing the 
results between two different frameworks. 

This paper is organized as follows. In Section [2] we ex- 
plain the shattering models and pick out some major param- 
eters. In Section [3] we examine the variation of the grain 
size distribution due to changes of the major parameters. In 
Section [4] we discuss our results and implications for the 
evolution of grain size distribution. In Section|5] we give our 
conclusions. 

2. Overview of the models for shattering 

As mentioned in Introduction, we use the framework in 



H10 to test the robustness of a shattering model. H10 con- 
sider the evolution of grain size distribution through shatter- 
ing by solving a shattering equation, which treats the grain- 
grain collision rate and the redistribution of shattered frag- 
ments in the grain size distribution. They assume that grains 
collide with each other under the velocity dispersions in- 
duced by dynamical coupling with interstellar turbulence. 
The production of fragments in grain-grain collisions is the 
key process that we focus on below. 

The dust grains are assumed to be spherical with radius a. 
Thus, the grain mass m is related to the grain radius by 

4 o 

m = -7TO p gri (1) 

where p gI is the material density of the grain. The grain size 
distribution, n(a), is defined so that n(a) da is the number 
density of grains with radii between a and a + da. 
2.1 Initial grain size distribution 

The size distribution of grains ejected from SNe II into the 
ISM is adopted from N07. This size distribution is used as 
the initial condition for the calculation of shattering. N07 
treated dust nucleation and growth in a SN II, taking into ac- 
count the dust destruction by kinetic and thermal sputtering 
in the shocked region. Thus, the grain size distribution cal- 
culated by N07 is regarded as that ejected from SNe II to 
the ISM (see also Bianchi and Schneider, 2007). Following 
H10, we adopt 20 M Q as a representative progenitor mass, 
and the unmixed case in which the original onion-like struc- 
ture of elements is preserved, since the extinction features of 
carbon and silicon, which are major grain components in the 
unmixed case, are consistent with observations (Hirashita et 
al., 2005; Kawara et al., 2010). The formed grain species are 
C, Si, Si0 2 , Fe, FeS, A1 2 3 , MgO, MgSi0 3 , and Mg 2 Si0 4 . 
According to N07, small grains with a < 0.02 /im are 
trapped in the shocked region and are efficiently destroyed 
by thermal sputtering if the ambient hydrogen number den- 
sity, nn, is larger than 0.1 cm~ 3 . In this paper, we adopt 
fin = 1 cm~ 3 (Section 12.21 ). As shown in N07 (see also 
H10), most of the grains have radii 0.1-1 /im. The initial 
grain size distributions are also shown later in Fig. |2] 

The normalization of the grain size distribution is deter- 
mined as follows. We consider a solar metallicity environ- 
ment as a representative case. (For other metallicities, the 
timescale of shattering just scales with the inverse of metal- 
licity.) We parameterize the dust abundance by the oxy- 
gen abundance. The oxygen mass produced by a SN II of 
20 M progenitor is mo = 1-58 M Q according to Umeda 
and Nomoto (2002), whose data were adopted by N07. The 
solar oxygen abundance is assumed to be Zo,© = 5.6 x 10~ 3 
by mass ratio (Lodders, 2003). Thus, for the solar oxygen 
abundance, the corresponding dust-to-gas ratio is given by 
Vq = Zq /mo > where is the dust mass ejected from 
a single SN II (0.14 M Q ; N07; H10). Here we assume that 
both oxygen and dust are supplied from SNe II. The grain 
size distribution is normalized so that the total grain mass 
density integrated for all the size range, pdust, is equal to 
lAnumaVo, where toh is the mass of hydrogen atom, and 
the factor 1.4 is the correction for the species other than hy- 
drogen. 
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Table 1. Summary of experimental properties. 



Species 


Pgr 


CO 


s 


Kshat 


Pi a 




(g cm" 3 ) 


(kms -1 ) 




(kmT 1 ) 


(dyn cm~ 2 ) 


Silicate 


3.3 


5 


1.2 


2.7 


3 x 10 11 


Carbonaceous dust 


2.2 


1.8 


1.9 


1.2 


4 x 10 10 



a Pi is varied in Models D and E in Table|2] 



Note: All the quantities are determined experimentally. See Jones et al. (1996) and references therein. 



As mentioned in Introduction, additional contribution 
from AGB stars does not change the following results signif- 
icantly as long as AGB stars also supply large (a > 0.1 /xm) 
grains (see also Asano et al., 2013). In other words, the use 
of N07's results for the initial condition is aimed at investi- 
gating the important role of shattering in the production of 
small grains from large grains. 
2.2 Shattering 

We calculate the shattering processes in the same way as 
in H10 by solving the shattering equation. The collision 
frequency between grains with various sizes is determined 
by grain-size-dependent velocity dispersions. The shattering 
equation used in H 10 is based on Hirashita and Yan (2009) 
(originally taken from Jones etal., 1994, 1996). Shattering is 
assumed to take place if the relative velocity between a pair 
of grains is larger than the shattering threshold velocities, 
Wshat (2.7 and 1.2 km s^ 1 for silicate and graphite, respec- 
tively; Jones et al., 1996). The results are not sensitive to 
w s hat as long as the grain velocities driven by turbulence is 
much larger than the threshold, but are rather sensitive to the 
hardness of the grain materials as shown later. 

We consider a range a m ; n = 3 x 10~ 8 cm (3 A) and 
ctmax = 3 x 10~ 4 cm (3 fim) for the grain radii. The 
grains passing through the boundary of the smallest radius 
are removed from the calculation. Most grain models adopt 
a minimum grain radius of a few x 10~ 8 cm (Weingartner 
& Draine, 2001; Guillet et al., 2009) although applying bulk 
material properties to such small grains could cause a large 
uncertainty. However, even if we adopt a m j n = 10~ 7 cm, the 
results does not change, except that the grain size distribution 
is truncated at 10 -7 cm. This is because shattering occurs in 
a top-down manner in the grain sizes and the production of 
grains with a ~ a few x 10~ 8 cm by shattering is never 
enough for these small grains to have a large contribution to 
the total shattering rate. 

Although nine grain species are predicted to form (Sec- 
tion |2Tj, the material properties needed for the calculation 
of shattering are not necessarily available for all the species. 
Thus, we divide the grains into two groups: one is car- 
bonaceous dust and the other is all the other species of dust 
(called "silicate"), and apply the relevant material quantities 
of graphite and silicate, respectively. The material properties 
of silicate and graphite are taken from Jones et al. (1996) and 
summarized in Table Q] 

Since the shattering equation is general enough, the major 
uncertainties can be produced by the treatment of shattering 
fragments. Thus, we explain how to treat shattering frag- 
ments below. 

We consider two colliding grains whose masses are mi 



and ra-i (the former grain is called target grain), and estimate 
the total fragment mass in ra\ as a result of this collision. 
Note that we consider a collision between m-i and m,\ again 
in the calculation and consider the fragments in mi (that is, 
we consider the same collision twice to treat the fragments of 
each colliding grain). The total fragment mass is determined 
by the mass shocked to the critical pressure in the target 
(Jones et al, 1996; Hirashita and Yan, 2009): 



M ej _ M _ f M (l + 2R) 1 / Ml 
m 2 /M m 2 2(1 + ft) 9 / 16 1/9 {^Mj 



8/9 



, (2) 



where M j is the total mass of fragments ejected from mi, 
M is the mass shocked to the critical pressure (Pi), above 
which the solid becomes plastic (that is, the deformation 
destroys the original lattice structures), Jm is the fraction 
of shocked mass that is eventually ejected as fragments (the 
rest remains in the grain), 1Z = 1 in the collision between the 
same species (we only consider collisions between the same 
species for simplicity), M r = v/cq (cq is the sound speed of 
the grain material), Mi is the Mach number corresponding 
to the critical pressure P\ : 



Mi = 



20i 



l + (l + 4s</>i) 1 /2 ! 



(3) 



where </>i = Pi/(/O gr c§) and s is a dimensionless material 
constant that determines the relation between the shocked 
velocity and the velocity of the shocked matter. For conve- 
nience, we define function a as 



a{M) 



o.mis + M- 1 -0.11) 



1.3 



M- 



1 



(4) 



and evaluate a\ = a(M\) and ay = er(.M r /(l + 7Z)) 
in Eq. ©. If M is larger than half of the grain mass, we 
assume that the whole grain is fragmented; i.e., M C j = mi. 
This case is called catastrophic disruption. Otherwise (i.e., 
for M < mi/2), only a fraction of mi (JmM) is ejected 
as fragments. This case is called cratering. The relation 
between the total volume (V) of shattered fragments and 
the critical pressure is roughly given by the total energy 
P\V = constant. Indeed, the above equations tell us that 
M oc V is inversely proportional to Pi, if 4s0i > 1, 
which is true for most of the cases considered in this paper 
12 ~ ^ so M ei oc 0^ 8/9 ex P~ 8/9 ). Since the 



'■'J 



{Ml ~ 

shattering efficiency is proportional to M e j, the timescale of 
shattering is regulated by Pi . Considering that the material 
properties of astronomical grains are uncertain, we vary Pi 
(Section l2A2l . 

The total fragment mass il/ C j is distributed with a grain 
size distribution oc a~ at . Jones et al. (1996) argue that 
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aif = 3.0-3.4, based on an analysis of the flow associated 
with the formation of a crater (see also Takagi et al., 1984; 
Nakamura and Fujiwara, 1991; Nakamura et al., 1994; Taka- 
sawa et al., 201 1 for experimental results). Unless the catas- 
trophic disruption occurs, we assume that the mass mi — M e j 
remains as a single dust grain and is distributed in an appro- 
priate bin in the numerical calculation. To avoid unnecessary 
complexity caused by the choice of the upper and lower radii 
of the fragments, we adopt simpler forms for the smallest and 
largest fragments than HI 1, following Guillet et al. (2009): 



a /max = (0.0204/) 1 / 3 a 1 



(5) 



where a\ is the radius of the grain m\ (i.e., the target grain), 
/ = Mej/mi (shattered fraction of the target grain). The 
minimum radius is assumed to be a/ >m in = a min = 3 x 
10~ 8 cm. If a/ imax < a m j n , we remove the fragments from 
the calculation. 

The relative velocity in the collision is estimated based 
on the grain velocity dispersion as a function of grain radius 
a. We used the calculation by Yan et al. (2004), who con- 
sider the grain acceleration by hydrodrag and gyroresonanc e 
in magnetohydrodynamic turbulence, and calculate the grain 
velocities achieved in various phases of ISM. Among the 
ISM phases, we focus on the warm ionized medium (WIM) 
to investigate the possibility of efficient shattering in actively 
star-forming environments. We adopt nn = 1 cm~ 3 for the 
hydrogen number density of the WIM. The resulting grain 
size distribution is not sensitive to «h (H10). Indeed, if 
riH is large, the grain-grain collision rate (i.e., the shatter- 
ing efficiency) rises, while dust grains are more destroyed in 
SNe (i.e., the dust-to-gas ratio is lower). Because of these 
compensating effects, the resulting grain size distribution is 
not sensitive to the gas density. We adopt gas temperature 
T = 8000 K, electron number density n c = nn, Alfven 
speed Va = 20 km s _1 and injection scale of the turbulence 
L = 100 pc. For grains with a > 0.1 /im, where most of 
the grain mass is contained in our cases, the grain velocity 
is governed by gyroresonance. Since there is still an uncer- 
tainty in the typical grain radius (a c ) above which gyrores- 
onance effectively works, we also address the variation of 
results by a c (Section[2A4|). 

The grain velocities given above are velocity dispersions. 
In order to estimate the relative velocity (V12) between the 
grains with mi and m^, whose velocity dispersions are v\ 
and t>2, respectively, and we average the size distributions of 
fragments for four cases V12 = V\ + 1*2, \vi — V2\, fi, and 
i'2 to take the directional variety into account. See H10 for 
more details. 

2.3 Duration of shattering 

Since we consider shattering in the WIM, it is reasonable 
to assume that the shattering duration (denoted as i) is deter- 
mined by a timescale on which the ionization is maintained. 
A typical lifetime of ionizing stars is < 10 Myr (Bressan et 
al., 1993; Inoue et al., 2000). We basically adopt t = 10 
Myr in this paper, but we also consider longer durations for a 
longer starburst duration or shattering over multiple star burst 
episodes. 

2.4 Parameters 

Based on the discussions in the previous subsection, we 
change the following parameters related to the fragments. 



The models are summarized in Table [2] where Model A 
adopts "fiducial" values for the parameters. 

2.4.1 Ejected fraction of shocked material How 

much fraction of the shocked material is finally ejected as 
fragments is uncertain. This factor is denoted as /m in 
Eq. (0. According to Jones et al. (1996), ]m is between 
0.3 and 0.6, while H10 adopted 0.4. Thus, we adopt 0.4 as 
a fiducial value (Model A) and also examine 0.3 (Model B) 
and 0.6 (Model C). 

2.4.2 Critical pressure Considering the uncertainty in 
the actual materials of astronomical grains, we vary P\ by 
an order of magnitude as Models D and E in Table [2] The 
standard value adopted from Jones et al. (1996) is P\ = 
3 x 10 11 and 4 x 10 10 dyn cm~ 2 for silicate and graphite, 
respectively (Table [T] Model A in Table |2). We adopt three 
times smaller and larger values for Pi in Models D and 
E, respectively. We call grains with low/high Pi soft/hard 
grains. 

There is another possible way of including the critical 
pressure in the model based on Kobayashi and Tanaka (2010, 
hereafter KT10). Jones et al. (1996) constructed the frag- 
mentation model based on shock wave propagation caused 
by a collision (see also Mizutani et ai, 1990). On the other 
hand, Holsapple (1987) provided a scaling formula for col- 
lisional outcomes based on laboratory experiments of his 
group (the review of the scaling formula is seen in Holsapple, 
1993). In this formula, M e j is proportional to 



IP : 



for (p <C 1, where 



1 mim 2 2 
E imp = - ■ v 

2 mi + mi 



(6) 



(7) 



is the impact energy between mi and mi, and is the 
specific impact energy at M e j = mi/2. Since M e j ~ mi 
for (p ^> 1, KT10 simply connected the cases for ip <C 1 and 
ip>l, providing a simple formula for M C j as 



M ei = - 
J 1 



9 



-m\. 



(8) 



Now we relate to P\. The dominant channel for 
the small grain production is cratering of large grains by 
small grains. If we assume that m\ ^> m^ (and as long 
as v 2 is not much larger than unity), we obtain M e j ~ 
m2V 2 /(2Q^ ) ). On the other hand, if we approximate Eq. 
© as M C j/m 2 ~ A(p gI v 2 / Pi), A ~ 1 for the range of 
quantities in this paper. Therefore, if the two models are 
equivalent, ~ Pi/(2p gr ). 

Thus, we examine KTlO's formulation as follows. We 
adopt Eq. © instead of Eq. (f2j) for the total mass of frag- 
ments with the same fragment size distribution (oc a~ Qf with 
the same upper and lower grain radii). The value of <p is eval- 
uated by Eqs. © and Q, where is given by 



Pi 

2p gI 



(9) 



This equation indicates that, if the energy per unit volume 
given to the grain exceeds ~ Pi/p gI , the major part of the 
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Table 2. Models. 



Name 


species 


fu 


-Pi 

(dyn cm -2 ) 


Of 


grain velocities 


A" 


silicate 


0.4 


3 x 10 11 


3.3 


Yan et al. (2004) 




graphite 


0.4 


A 1 nlO 

4 x 10 


3.3 


Yan et al. (2004) 


B 


silicate 


0.3 


3 x 10 11 


3.3 


Yan et al. (2004) 




graphite 


0.3 


4 x 10 


3.3 


Yan et al. (2004) 


C 


silicate 


0.6 


3 x 10 11 


3.3 


Yan et al. (2004) 




graphite 


0.6 


4 x 10 


3.3 


Yan et al. (2004) 


D" 


silicate 


0.4 


1 x 10 11 


3.3 


Yan et al. (2004) 




graphite 


0.4 


1.3 x 10 10 


3.3 


Yan et al. (2004) 


c a 
r. 


silicate 


n 1 


y x iu 




i an et al. \A\jV)<+) 




graphite 


0.4 


1.2 x 10 11 


3.3 


Yan et al. (2004) 


F 


silicate 


0.4 


3 x 10 11 


2.3 


Yan et al. (2004) 




graphite 


0.4 


4 x 10 10 


2.3 


Yan et al. (2004) 


G 


silicate 


0.4 


3 x 10 11 


4.3 


Yan et al. (2004) 




graphite 


0.4 


4 x 10 10 


4.3 


Yan et al. (2004) 


H 


silicate 


0.4 


3 x 10 11 


3.3 


Eq. (Toll 




graphite 


0.4 


4 x 10 10 


3.3 


Eq. (H 



"We also examine KTlO's formulation for fragments (Section l2.4.2t . 



grain is disrupted. We test Models A, D, and E by adopting 
KTlO's formulation. 

Note that we have adopted the approximation, Ajf C j ~ 
p gI v 2 /Pi. This is a good approximation for 4scf>i > 1 
((pi — Pi/(p gr CQ); see Section l2~2l . However, if <\>\ <C 
1, M cj /m 2 ~ (p gr w 2 /Pi)/0i for Eq. ©, so that Q* = 
(Pi/2p gr )0i gives a better approximation for <pi <C 1. In 
other words, if <fii is significantly smaller than 1 (i.e., for soft 
dust), Eq. (0 overestimates by a factor of </>i . Thus, if we 
calculate shattering by adopting Eq. grains are shattered 
less compared with Jones et al.'s case for soft dust. In this 
paper, because it is easy to interpret Eq. (|9), we simply adopt 
it for KTlO's formulation. 

2.4.3 Size distribution of shattered fragments The 
power index of the size distribution of shattered fragments, 
af, may be reflected in the grain size distribution after shat- 
tering. The canonical value adopted in H10 is at = 3.3, 
which is based on Jones et al. (1996). Here we also examine 
the cases of af = 2.3 and 4.3 (Models F and G, respectively), 
considering a wide range obtained in experiments (Takasawa 
et al, 2011). 

2.4.4 Grain velocities Grain motions are mainly 
driven by gyroresonance in the diffuse medium (Yan et al., 
2004). Fig. Q] presents the grain velocities as a function of 
grain radius in the WIM (Section 12.2b . Large grains with 
a > 10~ D cm are accelerated to ~ 20 km s^ 1 (i.e., simi- 
lar velocities to the Alfven velocity). Since larger grains are 
coupled with larger-scale turbulence, the motions of larger 
grains are less affected by the damping by dissipation, which 
occurs on small scales. The dissipation scale depends on the 
physical condition of the ISM, especially on the ionization 
degree. Thus, the critical grain size, above which the grains 
are efficiently accelerated by gyroresonance, varies by the 
physical condition of the ISM and the power spectrum of the 
turbulence. Because of such a variation of the critical grain 



size, it is worth investigating how the shattered grain size 
distribution changes as a result of the variation of the critical 
grain size. 

To make clear the dependence of the critical grain size 
denoted as a c , we also examine a model in which the size 
dependence of the grain velocity dispersion is simplified as 

v(a) = / %yroifa>a C) 
w [ ifa<a c . 

An example with a c — 10~ 6 cm is shown in Fig.Q] We fix 
Vgyro = 20 km s _1 . As long as Ug yro > w s h, the variation of 
Ugyro changes the timescale of shattering, which is roughly 
scaled with u~i„. We examine a c — 10~ 5 and 10~ 6 cm. 

gyro 

This treatment of velocities is labeled as Model H. 
3. Results 

3.1 Ejected fraction of shocked material 

We examine the effect of Jm on the grain size distribution. 
In Fig. |2] we show the grain size distributions at t = 10 Myr 
for Models A, B, and C. In this paper, the grain size distribu- 
tions are presented by multiplying a 4 to show the mass dis- 
tribution in each logarithmic bin of the grain radius. We ob- 
serve that the production of small grains is the most efficient 
for the largest Jm- This is simply because a larger fraction of 
the shocked material is ejected as fragments. However, the 
abundance of the smallest grains is not simply proportional 
to Jm, since the catastrophic disruption does not depend on 
/m- The maximum ratio of n(a) between Jm = 0.3 and 0.6 
is 1 .3 for silicate and 1 .6 for carbonaceous dust. 

3.2 Critical pressure 

We examine the dependence on the critical pressure (Pi). 
The adopted values for the fiducial case (Model A) is shown 
in Table|2l while smaller and larger values for Pi are also ex- 
amined (Models D and E, respectively). In Fig. [3] we com- 
pare Models A, D, and E at t = 10 Myr. We find that the ef- 
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Fig. 1. Grain velocities v as a function of grain radius a. Two grain species, (a) silicate and (b) carbonaceous dust, are shown. The solid lines indicate 
the velocity calculated by Yan et al. (2004) with the physical condition of the WIM described in the text. The dashed line shows the simplified velocity 
model (Eg. 1 101 with a c = 10 -6 cm. The dotted line shows the shattering threshold velocity. 




Fig. 2. Grain size distributions presented by multiplying a 4 to show the mass distribution in each logarithmic bin of the grain radius. The solid, dashed, 
and dot-dashed lines show the results at t = 10 Myr for Models A, B, and C (/m = 0.4, 0.3, and 0.6), respectively. The dotted line presents the 
initial grain size distribution before shattering. Two grain species, (a) silicate and (b) carbonaceous dust, are shown. We also show a slope of 0.5 [i.e., 
n(a) oc a -3 - 5 ]. 



feet of Pi is significant. Carbonaceous dust in Model D has 
the lowest Pi, so that the largest grains are efficiently shat- 
tered into smaller sizes. The mass loss rate of large grains 
with mass m is proportional to P x 8 ^ 9 (Section [2.2b ; hence 
a low Pi increases the number density of small (a < 10~ 6 
cm) grains. On the other hand, collisional cascades deter- 
mine the slope of the grain size distribution as —3.5, inde- 
pendent of Pi (Hellyer, 1970; Tanaka et al, 1996; KT10). 
The slope is therefore given by —3.5 for a < 10~ 6 cm, a 



value (—3.5) consistent with one derived observationally for 
the dust grains in the Milky Way and the Magellanic Clouds 
(Mathis et al, 1977; Pei, 1992). 

To show the dependence on the treatment of total crater- 
ing volume, we also show the results calculated by the sim- 
ple formulation of KT10 (see Section [2A2] for the formula- 
tion) in Fig. [3] We observe that KTlO's formulation is in fair 
agreement with our formulation based on Jones et al. (1996) 
for the fiducial (Model A) and hard dust (Model E), but the 
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Fig. 3. Same as Fig.[2]but for Models A, D, and E. The thick solid, dashed, and dot-dashed lines show the results at t = 10 Myr for Models A, D, and 
E (fiducial, hard, and soft grains), respectively, while the thin lines adopt the treatment of fragments by KT10 (the same line species shows the same 
model). The dotted line presents the initial grain size distribution before shattering. Two grain species, (a) silicate and (b) carbonaceous dust, are shown. 



discrepancy is relatively large for the soft dust (Model D) 
because Eq. (0 is not a good approximation for the soft dust 
(see Section 12.4.21 for detailed discussions). The maximum 
ratio between the two formulations in Model D is 1.8 for sil- 
icate and 2.2 for carbonaceous dust. For carbonaceous dust, 
the maximum difference occurs around the complicated fea- 
ture around a ~ 10~ 5 cm, but if we focus on the small grain 
production at a < 10~ 6 cm, the ratio is 1.4. Thus, the uncer- 
tainty caused by the different formulations is just comparable 
to that caused by the difference in fuj- If we adopt the stan- 
dard values for Pi (i.e., Model A), it is concluded that the 
small grain production by shattering is described by the sim- 
ple formulation by KT10 as precisely as more complicated 
framework of Jones et al. (1996). 
3.3 Size distribution of shattered fragments 

As we mentioned in Section 12.4.31 we assume that the 
size distribution of shattered fragments can be described by 
a power law function with power index ctf . In Fig. |U we 
show the results for different values of ctf (Models A, F, 
and G). We observe that the difference in ctf affects the 
resulting grain size distributions. However, we also see that 
the difference between ctf — 3.3 and 2.3 (Models A and F) is 
relatively small especially for silicate at a > 10~ 6 cm, which 
implies that the grain size distribution does not necessarily 
becomes n(a) cx a~ Qf . The grain size distribution in Model 
G at a > 10~ 6 cm is near to the initial grain size distribution, 
since the shattered grains are predominantly distributed at the 
smallest (a <C 10 -6 cm) sizes. 

In order to examine more clearly how the grain size dis- 
tribution evolves depending on acf, we show the results of 
Models F and G for t = 10, 20, 40, and 80 Myr in Fig. 
The dip around n~2x 10~ 7 cm, which is clearly seen at 



t = 80 Myr, is because of the velocity exceeding the shat- 
tering threshold around this grain radius (Fig. [T]. This non- 
monotonic behavior of the velocity is due to the complexity 
of the grain charge as a function of grain size (Yan et al., 
2004). 

There is a qualitative difference in the evolution of grain 
size distribution between at = 2.3 and 4.3. For at = 2.3, 
the slope of the grain size distribution at a < 10~ 6 cm 
changes, while for ctf = 4.3, it evolves with a constant 
slope [n(a) cx a~ at ] at small sizes. In the former case, 
the slope of the grain size distribution seems to approach 
n(a) cx a~ 3 5 . The variation of slope occurs in a top-down 
manner; for example, for carbonaceous dust, the grain size 
distribution is nearly n(a) cx a -3-5 between a = 10~ 6 
and 10 -5 cm at 40 Myr, while it approaches n(a) cx a -3 ' 5 
at smaller sizes later at 80 Myr. This top-down behavior 
is because of the nature of shattering, which disrupts large 
grains into a lot of small pieces. This trend is also seen for 
ocf = 3.3. Generally, if ctf < 3.5, the slope of grain size 
distribution tends to approach the MRN value (—3.5). On the 
contrary, if ctf > 4 (see also Section l4~2l . the resulting grain 
size distribution is determined by the grain size distribution 
of shattered fragments. Indeed, if ctf > 4, the smallest 
particles, which are too small to disrupt large grains around 
the peak (i.e., a > 10~ 5 cm), have the dominant contribution 
to the total fragment mass. Moreover, we remove the grains 
reaching the smallest-size bin (a = 3 x 10~ 8 cm), and these 
grains do not affect the subsequent evolution of grain size 
distribution. Thus, for ai > 4, the fragments just keep 
their grain size distribution and eventually lost in reaching 
the smallest grain size. 

We note that continuous shattering for such a long time as 
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Fig. 4. Same as Fig.[2]but for Models A, F, and G, in order to investigate the effect of cif (power index of the grain size distribution of shattered fragments). 
The solid, dashed, and dot-dashed lines show the results at t = 10 Myr for Models A, D, and E (of = 3.3, 2.3, and 4.3), respectively. The dotted line 
presents the initial grain size distribution before shattering. 




Fig. 5. Time evolution of grain size distribution for Models F and G in Panels (a) and (b), respectively, for carbonaceous dust. The thin solid, dashed, 
dot-dashed, and thick solid lines represent t = 10, 20, 40, and 80 Myr, respectively. The dotted line presents the initial grain size distribution before 
shattering. 



80 Myr is not applicable to the real WIM, since the lifetime 
of the WIM is probably shorter than 80 Myr. Yet, the above 
calculations for the long shattering durations show that, if 
grains are repeatedly shattered with ai < 3.5, the final grain 
size distribution can approach a power law with an index of 
~ —3.5. Therefore, shattering is a promising mechanism to 



produce an MRN-like grain size distribution. This issue is 
further discussed in Section PPl 
3.4 Grain velocities 

To understand the effect of grain velocities, we examine 
Model H, in which we use Eq. ( flOb and vary a c . In Fig. [6] 
we show the grain size distributions at t = 10 Myr for 
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Fig. 6. Same as Fig.[2]but for Model H with various a c . The solid, and dashed lines show the results for a c = 10 5 and 10 6 cm, respectively. The 
dotted line presents the initial grain size distribution before shattering, and the dot-dashed line shows the results for Model A as a reference. 



a c = 10~ 5 and 10~ 6 cm. We observe that the difference 
is clear at small grain sizes: because small grains are also 
shattered efficiently for small a c , the abundance of smallest- 
sized grains is larger for smaller a c . Since carbonaceous dust 
is softer than silicate, the effect of a c has a larger imprint 
on the grain size distribution of carbonaceous dust. Indeed, 
we clearly see the effect of a c in the grain size distribution 
of carbonaceous dust in Fig. [6] (right panel). Thus, if small 
(a < 10~ 6 cm) grains also acquire a larger velocity than 
the shattering threshold, it is predicted that the grain size 
distribution may not be approximated by a simple power law. 

In Fig. |6l we also show the results of Model A as a refer- 
ence. Since the grain velocities at a < 10~ 5 cm drops below 
20 km s- 1 for Model A (Fig. ©, Model H with a c = 1CT 5 
cm produce a more similar grain size distribution to Model 
A than Model H with a c = 10~ 6 cm. Thus, the effect of a c 
(i.e., the minimum size of grain size accelerated to the maxi- 
mum velocity) is more important than the detailed functional 
form of v(a) at a < a c . 

4. Discussion 
4.1 Timescales 

In order to interpret the results above, it is useful to es- 
timate a typical timescale of grain-grain collision. A grain 
with a = ai collide with a grain whose typical radius is a 2 
with a frequency of 



1 



7"coll 



7r(ai + a 2 ) 2 n(a 2 ) Aa 2 vi 2 ' 



(11) 



where Aa 2 is the radius interval of interest. Note that 
n(a 2 ) Aa 2 has a dimension of number density (Section 
[2}. By using logarithmic size interval A log a 2 , Aa 2 ~ 
a 2 Aloga 2 . If we consider a logarithmic size bin so that 



A log a 2 is of order unity, we can replace Aa 2 with a 2 . We 
estimate the typical timescale of small grain production by 
adopting the timescale on which a large grain is hit by a small 
grain (ai 3> a 2 ). In such a case, we obtain 



Tcoll 



tfn(a 2 )a 2 ^i 



[n(a 2 )a 



I]— v i 
a 2 



5.0 



"2 



io- 



20 km s~ 



u cm 
-l 

Gyr 



n(a 2 )au 



-i -l 



10 



-30 



(12) 



Thus, we find that the timescale of small grain production 
becomes shorter if the typical size of colliding grains (a 2 ) 
becomes small and the abundance of such grain [n(a 2 )] in- 
creases. For the initial grain size distribution, the grain 
radius is concentrated in a relatively narrow range around 
~ 3 x 10 -5 cm (Section 12.11 Fig. [2] and subsequent fig- 
ures also show the initial grain size distribution), so that 
a 2 /ai ~ 1/10 and a 2 ~ 3 x 10 -6 cm. Under such a con- 
dition, T co n ~ 1.5 x 10 8 yr, which is long compared with 
the typical duration treated in this paper (10 Myr). Thus, the 
large grains are hardly affected by shattering which is con- 
sistent with the fact that the grain size distribution around 
10~ 5 cm is intact at 10 Myr (Fig. |2j. However, after a small 
amount of shattering, a lot of small grains are produced, and 
if grains with a < 10~ 6 cm achieve n(a)a 4 ~ 10~ 30 , the 
collision timescale with such small grains for a large grain 
(ai ~3x 10~ 5 cm) becomes an order of 10 Myr. 

The timescale on which a large grain is lost because of 
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repeated collisions of small grains can be estimated by 

Tdisr ~ (mi/M oj )r coU . (13) 

This timescale is called disruption timescale. Here we show 
that the timescale of the loss of large (a ~ 3 x 10~ 5 cm) 
grains is consistent with the disruption timescale. For an- 
alytical convenience, we adopt KTlO's formulation, which 
is roughly equivalent with Jones et al. (1996)'s formula- 
tion through Eq. (|9). Since we are interested in the disrup- 
tion of large grains by collisions with small grains, mi 3> 
7Ti 2 holds. Thus, from Eqs. ® and (0, we obtain ip ~ 
(012/ 'oi) 3 '(pgrV 2 1 Pi). If we adopt a\ ~ 3 x 10~ 5 cm, 
a2 ~ 3 x 10~ 6 cm, and v ~ 20 km s -1 , we obtain y> ~ 0.22 
for carbonaceous grains. Thus, M c -Jmi ~ 0.18 from Eq. 
®, and we finally obtain Tdisr ~ 5.5t co h. Recalling that 
T co n is of an order of 10 Myr, the disruption timescale for the 
large grains is a few tens of Myr. This timescale is consistent 
with the loss of the grains around a~3x 10~ 5 cm in a few 
tens of Myr. 

4.2 Sensitivity to the parameters 

The resulting grain size distribution is sensitive to various 
parameters in different ways. As expected, shattering is more 
efficient for a smaller value of Pi simply because the grains 
are softer. Since the cratering volume is approximately pro- 
portional to 1/Pi (Jones et al., 1996), it is expected that 
similar grain size distributions are realized if we keep Pi jt 
constant. To examine if this is the case, we compare the re- 
sulting grain size distributions for constant P\/t in Fig. [7] 
We examine (Pi, t) = (3 X 10 11 dyn cm" 1 , 10 Myr), (1 x 
10 11 dyn cm -1 , 3.3 Myr) (i.e., soft dust and short shattering 
duration), and (9 x 10 11 dyn cm -1 , 30 Myr) (i.e., hard dust 
and long shattering duration) for silicate, and (Pi , t) = (4 X 
lO^dyncm -1 , 10 Myr), (1.3x 10 10 dyn cm -1 , 3.3 Myr), 
and (1.2 x 10 11 dyn cm -1 , 30 Myr) for carbonaceous dust. 
We observe that similar grain size distributions are indeed 
obtained within a factor of 3. The difference arises because 
the scaling of the cratering volume (or M e j) with P^ 1 is not 
perfect (Section l2~2l i. 

The power index of fragment size distribution, ct{, is 
simply reflected in the resulting grain size distribution for 
off = 4.3. However, if ai is smaller than 3.5 (e.g., the 
case with af = 2.3), the slope of grain size distribution at 
a < 10~ 5 cm continuously "steepens" and approaches to the 
slope consistent with MRN [n(a) oc a -3 5 ] (Fig. [4]). Accord- 
ing to the analytical studies of shattering by Hellyer (1970), 
the power index of grain size distribution after shattering has 
a steady state value of ~ —3.5 if the size distribution of shat- 
tered fragments is assumed to be a power law with an index 
smaller than 4 (i.e., a.{ < 4) (see also Tanaka et al., 1996), 
although we should note that the relative velocity between 
grains is a complicated function of grain sizes in our case. 
For af > 4, since the mass occupied by the smallest frag- 
ments is dominating, the size distribution after shattering is 
just dominated by the smallest fragments, keeping a power 
law index of — ot{ at small sizes. 

In general, larger grains have larger velocity dispersions if 
they are accelerated by turbulence, primarily because they 
are coupled with larger-scale turbulent motion which has 
larger velocity dispersions. Thus, even if the initial grain size 
distribution is biased to large grains with a > 10~ 5 cm, these 



grains are efficiently shattered to produce a large abundance 
of small grains. If these shattered grains have also larger 
velocities than shattering threshold, they collide with each 
other, accelerating the production of small grains. Thus, as 
shown in Fig. [6] the effect of a c (or the critical grain radius 
above which the grains have larger velocities than the shatter- 
ing threshold) can be seen in the grain size distribution, caus- 
ing a deviation from a monotonic power-law-like functional 
form. The imprint of a c is more clear in carbonaceous dust 
than in silicate dust, since the former species is softer. In- 
terestingly, Weingartner and Draine (2001) have shown that 
grain size distributions derived from the Milky Way extinc- 
tion curves are complicated for carbonaceous dust, which 
may be explained by imprints of a c . 

4.3 Implication for the evolution of grain size distribu- 
tion in galaxies 

Inoue (2011) shows theoretically that the overall dust 
mass in metal-enriched systems such as the Milky Way is 
regulated by the balance between dust growth in molecu- 
lar clouds and dust destruction by sputtering in SN shocks. 
While both these two mechanisms deplete small grains (Hi- 
rashita and Nozawa, 2012), shattering is a unique mecha- 
nism that reproduces the small grains. Indeed, Hirashita and 
Nozawa (2012) show that, unless we consider shattering, we 
underproduce the small grain abundance in the Milky Way 
(see also Asano et al., 2013). As shown in Section [3~3l shat- 
tering not only produces small grains but also has an inherent 
mechanism of steepening the grain size distribution; that is, 
even if the grain size distribution of shattered fragments has 
af smaller than 3.5, the resulting grain size distribution has 
a steeper dependence, approaching a slope consistent with 
the one derived by MRN [i.e., n(a) oc a -3 ' 5 ]. Thus, the 
robust prediction is that, if shattering is the dominant mech- 
anism of small grain production, the grain size distribution 
approaches an MRN-like grain size distribution as long as 
Off < 3.5. And if shattering is the dominant mechanism of 
small grain production, af > 4 is rejected, because the grain 
size distribution in the Milky Way ISM has a power index of 
-3.5 (MRN). 

In this paper, we have not included dust production by 
AGB stars. Even if the dust supply from AGB stars has a 
significant contribution, the results in this paper are not al- 
tered as long as the dust grains formed by AGB stars are 
biased to large (> 0.1 /im) sizes. Rather, the importance 
of shattering is emphasized because the additional dust pro- 
duction by AGB stars enhances the dust abundance, mak- 
ing grain-grain collision more frequent. Production of large 
grains from AGB stars is indicated observationally (Groe- 
newegen, 1997; Gauger et al., 1998; Norris et al, 2012) and 
theoretically (Hofner, 2008; Mattsson & Hofner, 2011). We 
also neglected grain growth in molecular clouds, which is 
suggested to be a major dust formation mechanism even at 
high redshift (Mattsson, 2011; Valiante et al, 2011). Grain 
growth, however, cannot be a supplying mechanism of small 
grains. 

For the shattering duration, since it is degenerate with the 
critical pressure of grains (hardness of grains), it is difficult 
to constrain it directly from the models. If grains experience 
shattering for several tens of Myr, the grain size distribu- 
tion possibly approaches a power law whose power index is 
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Fig. 7. Same as Fig.|3]but for constant Pi/t. The solid, dashed, and dot-dashed lines show the results for (a) (Pi, t) = (3 X 10 11 dyn cm , 10 Myr), 
(1 X 10 11 dyn cm" 1 , 3.3 Myr), and (9 X 10 11 dyn cm" 1 , 30 Myr), respectively, for silicate, and (b) (Pi, t) = (4 X 10 10 dyn cm" 1 , 10 Myr), 
(1.3 X 10 10 dyn cm" 1 , 3.3 Myr), and (1.2 X 10 11 dyn cm -1 , 30 Myr), respectively, for carbonaceous dust. The dotted line presents the initial 
grain size distribution before shattering. 



consistent with the MRN, if at < 3.5 (Fig. |5}. Since the 
timescale of ISM phase exchange is comparable or shorter 
than that required for the grain size distribution to approach 
the MRN size distribution (O'Donnell and Mathis, 1997), 
shattering may occur intermittently. Even in this case, it 
is expected that grains are relaxed into an MRN-like sim- 
ple power law if the total shattering duration reaches several 
tens of Myr. Such a duration for shattering is probable since 
it is much shorter than the grain lifetime (~ a few x 10 8 Myr; 
Jones etal, 1996). 

5. Conclusion 

Shattering is a viable mechanism of small grain produc- 
tion in the ISM both in nearby galaxies and in high-redshift 
galaxies. We have examined if grain size distributions pre- 
dicted by shattering models are robust against the change 
of various parameters and formulations. Because of the un- 
certainty in Jm (the fraction of the shocked material that is 
eventually ejected as fragments), the predicted grain size dis- 
tribution after shattering is uncertain by a factor of 1.3 (1 .6) 
for silicate (carbonaceous dust). We have identified Pi (the 
critical pressure above which the original lattice structure is 
destroyed) as the most important quantity in determining the 
timescale of small grain production, and confirmed that the 
same Pi jt (t is the duration of shattering) gives roughly the 
same grain size distributions after shattering within a factor 
of 3. 

A simpler and more intuitive formulation by KT10 is in 
good agreement with our model based on Jones et al. (1996) 
within a factor of 1.8 (1.4) for silicate (carbonaceous dust) 
if we focus on the small grain production at a < 10 ~ 6 cm. 
This is as small as the uncertainty caused by /m- Thus, as 



long as we have an uncertainty in /m, it is sufficient to adopt 
the simpler formulation by KT10. The size distribution of 
shattered fragments have minor effects as long as at < 3.5, 
since the grain size distribution is continuously steepened by 
shattering and become consistent with the MRN grain size 
distribution regardless of the value of at (< 3.5). 

The effect of the grain velocities as a function of grain 
radius can be seen more clearly in carbonaceous grains than 
in silicate, since the former species is shattered more easily. 
Thus, it is predicted that carbonaceous species has more 
complicated grain size distribution than silicate; in other 
words, the grain size distribution of carbonaceous dust can 
show some imprints of the grain velocity as a function of 
grain radius. 
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